Femtosecond dynamics of the collinear-to-spiral antiferromagnetic 

phase transition in CuO 

S. L. Johnson, 1 R. A. de Souza, 1 U. Staub, 1 P. Beaud, 1 E. Mohr-Vorobeva, 1 
G. Ingold, 1 A. Caviezel, 1 V. Scagnoli, 1 W. F. Schlotter, 2 J. J. Turner, 2 0. 
Krupin, 3 ' 2 W.-S. Lee, 4 - 5 Y.-D. Chuang, 6 L. Patthey, 1 R. G. Moore, 4 D. 
Lu, 5 M. Yi, 4 P. S. Kirchmann, 4 M. Trigo, 7 P. Denes, 6 Dionisio Doering, 6 
Z. Hussain, 6 Z.-X. Shen, 4 D. Prabhakaran, 8 and A. T. Boothroyd 8 

1 Swiss Light Source, Paul Scherrer Institut, 
CH-5232 VilUgen PSI, Switzerland 
2 The Linac Coherent Light Source, 
SLAC National Accelerator Laboratory, Menlo Park, CA 94025 USA 
^European XFEL GmbH Albert Einstein Ring 19, 22 607 Hamburg, Germany 
4 SIMES, SLAC National Accelerator Laboratory and 
Stanford University, Menlo Park, CA 94305, USA 
5 SSRL, SLAC National Accelerator Laboratory, Menlo Park, CA 94305, USA 

6 'Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
7 PULSE, SLAC National Accelerator Laboratory, Menlo Park, CA 94025 USA 
^Department of Physics, University of Oxford, 
Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, UK 



1 



Abstract 

We report on the ultrafast dynamics of magnetic order in a single crystal of CuO at a temperature 
of 207 K in response to strong optical excitation using femtosecond resonant x-ray diffraction. In 
the experiment, a femtosecond laser pulse induces a sudden, nonequilibrium increase in magnetic 
disorder. After a short delay ranging from 400 fs to 2 ps, we observe changes in the relative intensity 
of the magnetic ordering diffraction peaks that indicate a shift from a collinear commensurate 
phase to a spiral incommensurate phase. These results indicate that the ultimate speed for this 
antiferromagnetic re-orientation transition in CuO is limited by the long-wavelength magnetic 
excitation connecting the two phases. 

PACS numbers: 75.25.-j,75.30.Kz,75.78.Jp 
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The ability to manipulate magnetic order has become vital to modern data storage and 
processing technology. As this technology is pushed to handle ever increasing data rates, 
fundamental questions arise on the ultimate time scales required for changes in magnetic 
order and how it is possible to induce them in a controlled way p~H6]. In several ferromagnetic 
materials, experiments have demonstrated demagnetization on timescales of several hundred 
femtoseconds [7HTT]. While this addresses one aspect of the problem, another potential 
focus is the dynamics of phase transitions between different types of magnetic order [121 - 
[T5] . Of particular interest here are phase transitions involving antiferromagnetic spiral spin 
structures resulting from magnetic frustration, which often lead to multiferroicity [IS]. CuO, 
well-known as a starting material for high-T c cuprate superconductors, shows one important 
example of such a phase transition with some unusual characteristics. 

The atomic structure of CuO can be described as connected networks of zigzag Cu-0 
chains oriented along the [10—1] and [1 1] directions [T7j. Depending on the temperature, 
the spins on the Cu atoms take on different ordering patterns [T8TI20"] . Figure 1 shows the 
structure of the magnetically ordered phases of CuO, projected into the (0 10) plane. The 
low temperature phase has a collinear antiferromagnetic structure. It has a commensurate 
(CM) ordering wave-vector (0.5 -0.5) that forms a magnetic unit cell twice the dimensions 
of the structural unit cell in both the a and c directions. A multiferroic phase appears in the 
temperature range 213—230 K [2T]. In this phase the magnetic order is qualitatively different: 
it is a non-collinear antiferromagnet with an incommensurate (ICM) ordering wave-vector 
(0.506 —0.483) that is not pinned to the underlying lattice periodicity. The spiral spin 
structure of this phase results in a small ferroelectric moment. Interestingly, the temperature 
at which this induced multiferroicity exists in CuO is much higher than has been observed in 
other materials, a consequence of the unusually large magnetic superexchange interactions 
in copper oxides that may also have a connection to the high-T c superconductivity often 
observed in these compounds [21]. At temperatures above 230 K CuO loses long-range 3-D 
magnetic order. 

These changes in magnetic order are the result of frustration between ferromagnetic and 
antiferromagnetic exchange interactions among the spins on the Cu atoms. The strongest 
coupling in the system is a nearest-neighbor antiferromagnetic (AF) interaction along the 
[10-1] direction, with exchange interaction Jaf = 74 meV [22[ 123]. The large strength of 
this interaction results in 1-D collinear AF chains along [10 —1] in both the CM and ICM 
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phases. To explain the CM/ICM transition, Yablonskii [23] proposed a model that considers 
the energy E from effective nearest-neighbor and next-nearest-neighbor interactions among 
these AF chains: 



where S n refers to the spin of a Cu atom along the c axis with index n as shown in figure 1, 
Ji < is a nearest-neighbor ferromagnetic exchange interaction, J2 > is a next-nearest- 
neighbor antiferromagnetic exchange interaction, and the biquadratic exchange interaction 
/ > is assumed to be small relative to the others. The biquadratic interaction term acts 
to stabilize the CM phase at low temperature. By invoking a mean-field approximation and 
taking S as the spatially averaged value of the spin magnitude, the CM phase has a lower 
energy for S 2 > Jf/8IJ2- Since S decreases with temperature, at some point the value of 
S 2 drops below this threshold and the ICM phase becomes more stable. 

Although this particular role of a biquadratic energy term as the mechanism for the low-T 
stability of the CM phase may be unique to CuO, the role of thermal disorder in the phase 
transition is more general. In a mean- field treatment, temperature is essentially modeled as 
a change in the average value of an order parameter due to disorder. In multiferroic and 
other materials with competing interactions and order parameters this can result in order- 
order phase transitions. CuO presents one particular example of how this disorder can lead 
to a qualitative change in the nature of magnetic order. 

Temperature is only one way to control the level of spin disorder and thus the value of S. 
Alternatively, optical excitation of CuO above the charge transfer bandgap at 1.4 eV [25] 
can directly induce non-thermal spin disorder by exciting electron transfer. This allows us 
to use an intense, short pulse of light to nearly instantaneously drive the magnetic order 
toward the ICM phase. We can then explore the limiting time scales of the transition itself. 

We use femtosecond time-resolved resonant x-ray scattering to observe the effect of optical 
excitation on the CM and ICM magnetic structures. For x-rays with an energy tuned 
near the Cu L 3 edge at 930 eV, the scattering from magnetic contributions is strongly 
enhanced. The diffracted intensity at the magnetic ordering wavevector has no measurable 
structural component, but is instead a direct measure of the ordered component of the 
dipole magnetic moment and an orbital flux (toroidal moment) which coexist within the 
CuO plaquette [261 121] ■ The orbital flux component is an important ingredient in the 
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description of the pseudo-gap phase order parameter in cuprate superconductors [28]. By 
using < 10 fs pulses of x-rays in combination with synchronized 40 fs pulses of 800 nm light, 
we investigate the magnetic dynamics using a pump-probe method. 

The sample is a single crystal of CuO grown from high-purity CuO powder using the 
floating- zone technique [29]. Cylindrical rods of length 10 cm and diameter 12 mm were 
sintered for a three day period in a vertical furnace at 950° under oxygen. The crystal was 
then grown in an optical floating-zone furnace using a mixed oxygen/argon atmosphere. The 
surface was miscut at an angle of 10° with respect to the (0.5 -0.5) planes. 

The pump-probe resonant x-ray diffraction measurements were performed at the Linac 
Coherent Light Source (LCLS) [30J, using the soft x-ray materials science (SXR) instrument. 
The crystal was cooled to a temperature of 207 K in vacuum using a He-flow cryostat. To 
excite the sample, a 40 fs full-width-at-half maximum pulse of 800 nm wavelength light 
hit the sample with an horizontal angle of 64.1° and 7r-polarization with respect to the 
(0.5 -0.5) planes. Sub-10 fs pulses of x-rays probed the magnetic ordering. For these 
measurements the photon energy was tuned to optimize output at 930 eV. A variable-line- 
spaced plane-grating monochromator set this as the x-ray energy to a bandwidth of 1.2 eV. 
A set of Kirkpatrick-Baez mirrors focused the x-ray beam to a 300 micron diameter spot on 
the sample. The incident x-ray beam was collinear with the excitation beam and 7r-polarized 
with respect to the (0.5 -0.5) planes. A continuously moving optical delay line in the path 
of the pump beam set the relative delay of the two synchronized pulses. For each pulse, this 
timing was corrected by electron bunch time arrival monitors. 

A custom-built CCD camera detected the diffracted x-rays. The camera is able to isolate 
individual x-ray pulses from the 60 Hz repetition rate of the LCLS. The sample angle was 
tuned intermediate between the optimal Bragg conditions for the CM and ICM phases, which 
are separated by less than 1.5°. This allowed simultaneous measurement of both peaks. The 
peak positions did not change as a function of time delay, indicating that changes in the 
normal component of the reciprocal lattice vectors (due to, for example, lattice expansion) 
do not contribute to the measured data. To correct for intensity and energy fluctuations of 
the source, the overall image intensities were normalized by the drain current signal from a 
thin Al-backed foil placed in the x-ray beam before the sample. The intensity of the peaks 
for a given image was estimated by summing the image data along the short axis of the 
CCD image and then fitting the resulting curves to the sum of two Gaussians. The integral 
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of each Gaussian was then taken as the intensity. 

Figure 2 shows images of the magnetic diffraction signal at times before and after exci- 
tation. The two peaks show contributions from each of the two types of magnetic ordering. 
They coexist in different domains at this temperature due to the first-order nature of the 
phase transition. Because the estimated probe depth of 50 nm is less than the correlation 
length (approximately 100-200 nm), the measured domains are distributed predominantly 
along the surface of the sample. The majority of these domains are initially in the CM 
phase. After excitation, the intensities of both peaks decrease but the relative intensity 
shifts in favor of the ICM phase, indicating that the optically-induced spin disorder has 
indeed changed the average magnetic order. 

Figure 3a shows the time-dependence of the intensities for a ffuence of 28 mJ/cm 2 . When 
the optical pulse hits the sample, both peaks show a sudden decrease over a time scale of 
300 fs. The relative change in intensity here is the same for both peaks; we also observe 
this for other excitation ffuences. This sudden drop has a simple interpretation: it is a 
direct measure of the additional spin disorder in the system induced by optical excitation. 
This is analogous to the drop in intensity observed in structural diffraction peaks due to 
non-thermal lattice disorder [3TJ E2] • The fast speed indicates that the disordering occurs 
on a time scale faster than the estimated experimental time resolution of 300 fs. The exact 
mechanism of the laser-induced spin disorder is not unambiguously determined from the 
experimental results. 

At a time slightly later than the initial drop, the time dependence of the two peaks 
begins to differ. Figure 3b shows the time-dependence of the ratio between the ICM and 
CM intensities. These data show the relative contributions of each phase as a function of 
time after the onset of spin disorder. The curve superimposed on the data is a fit to a 
bi-exponential curve that begins at a delayed time t p . As shown in the inset to figure 3b, 
the fitted value of t p is significantly larger than zero for all measured excitation fluences. 
The delay is 2 ps for low excitation levels and becomes shorter with increasing fluence 
until reaching a "saturation" level of approximately 400 fs. Note that since t p is derived 
from a comparison of mathematically independent data sets, it is not limited by the time 
resolution of the experiment. After t p , the relative intensity of the incommensurate peak 
grows. Initially there is a fast component with a time constant of 1.1 ± 0.3 ps for the 
lowest fluence of 7 mJ/cm 2 , increasing slightly with fluence to a value of 1.86 ± 0.14 ps at 
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39 mJ/cm 2 . There is also a much slower time scale component with a time constant ranging 
from approximately 35-60 ps, becoming faster with increasing ffuence. In this context it 
is important to note that these fits considered only times up to 20 ps after excitation; the 
later time behavior is not quantitatively well described by this model. The ratio at 100 ps 
increases with excitation fluence [33]. At the maximum fluence of 39 mJ/cm 2 we observed 
an increase of the ICM/CM ratio by a factor of 2.8 over the initial ratio. 

In principle, there are at least two possible explanations for the difference in the time 
response between the two phases. One possibility is that the relaxation processes for the two 
magnetic phases differ intrinsically. The dynamics can arise from both changes in the orbital 
flux and the magnetic moment contributions to the magnetic scattering. Although orbital 
current is expected to couple more strongly to the lattice since it is related to inversion 
symmetry breaking, we expect that the much stronger interactions between magnetic dipole 
moments control the dynamics in CuO. The magnetic dipole disorder induced by optical 
excitation of the charge-transfer band is expected to create predominantly local exchange of 
spins among Cu atoms, resulting in a broad spectrum of magnetic excitations throughout the 
Brillouin zone. Interactions of these magnetic excitations both with each other and with the 
lattice and charge degrees of freedom will at some point result in a thermalized distribution 
of magnetic excitations in equilibrium with other energy subsystems in the crystal. Since 
the initial steps of this process involve relaxation of high-frequency magnetic excitations 
{hv ~ 150 meV) [23], it is not offhand clear why there should be a difference between 
two phases that mostly differ in their low-frequency, long-wavelength magnetic behavior. 
Measurements (see |33j) of the time-dependent intensity of the CM peak at 150 K where no 
ICM peak was observed show behavior that is qualitatively different: the intensity of the 
CM peak increases over a time scale of 2-3 ps. It thus appears that the decrease in the CM 
peak intensity observed closer to the phase transition temperature is strongly related to the 
transition itself. 

Alternatively, the difference in the time response can be the result of a shift in population 
of the CM to the ICM phase, possibly starting at the boundaries between pre-existing 
domains. The excitation-dependent delay may be explained qualitatively in terms of the 
configuration schematics in figure 4. In general, any change in structure can be described 
as motion along a set of coordinates associated with normal excitation modes. Although in 
fact the set of magnetic coordinates that connect the CM and ICM phases involves several 
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different modes, for the sake of simplicity we consider only changes along one coordinate 
that physically corresponds to an acoustic-branch magnetic excitation with wavevector q = 
1\CM ~~ ^CM ~ 0- This is a useful approximation since this is likely to be among the lowest 
frequency modes that mediate the change in magnetic order. Initially, the energy of the CM 
phase is lower than that of the ICM phase. The energy barrier between the CM and ICM 
phases may be roughly estimated as the magnitude of the biquadratic term in equation [T], 
evaluated as IS 4 in the mean field approximation for the CM phase. At low excitation 
fluences, the relative energies of the two phases shifts in favor of the ICM phase. The energy 
barrier is, however, still sizable. The system requires a non-zero amount of time to hop over 
the energy barrier and enter into the new phase. As the excitation fluence increases, the 
size of the energy barrier decreases until it reaches a level where it no longer appreciably 
impedes the speed of the transition. At these very high excitation levels, we might expect 
that the time scale is limited purely by dynamics: the time required for a coherent, 1/4 
oscillation of a low momentum q w spin wave. Inelastic neutron scattering measurements 
in the CM phase have estimated a period of 1.6 ps for such excitations [19] . Taking 1/4 
of this period yields 400 fs, in agreement with the delay observed at the highest excitation 
fluences. This is analogous to the situation with structural phase transitions, where time 
scales of a fraction of a phonon period are frequently cited as minimum time scales [34], 

After t p , the fast increase of the ICM/CM ratio can be considered as the first step of the 
actual phase transition. The subsequent increase of the ratio over several tens of picoseconds 
occurs at approximately time as a slow decrease in the overall diffracted intensity. We 
tentatively interpret this as an intermediate step in thermalization that further suppresses 
the spatially averaged S 2 and pushes the system even more toward the ICM phase. 

The concept of using intense optical pulses to control magnetic order has the potential 
for a wide variety of applications. We have shown here one possible way to accomplish a 
partial, transient change in average magnetic order from one type to another by using a 
40 fs duration optical pulse to initiate the process. In addition, we have identified the 
limiting time scales for how fast the first steps of such a change may occur. The observed 
changes in magnetic order may be determined by the growth of pre-existing domains; further 
studies with improved spatial resolution may help to clarify this. 

This research was carried out on the SXR Instrument at the LCLS, a division of SLAC and 
an Office of Science user facility operated by Stanford University for the U.S. Department 
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FIG. 1. Magnetic structures in CuO, projected onto the (0 1 0) plane [T71|35j. Left: the magnetic 
structure of the CM phase; right: that of the ICM phase. Red arrows indicate the orientation of 
the magnetic dipoles. 
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FIG. 2. Images of resonant x-ray diffraction. Upper panel: 1 ps before excitation. Lower panel: 
100 ps after excitation with a femtosecond laser pulse at an incident fluence of 39 mJ/cm 2 . The 
false-color scale for each image is normalized to better show the relative intensities. 
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FIG. 3. (a) Time-dependent normalized diffracted intensity for an excitation fluence of 28 mJ/cm 2 . 
Inset: longer time-scale behavior, (b) Relative phase population. The curve is a fit described in 
the text. Inset: dependence of the onset time t p on the excitation fluence. The curve is a guide to 
the eye. The dashed line shows the time for a 1/4 period of a long- wavelength magnetic excitation 
in the CM phase (400 fs). See [33] for additional data. 
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FIG. 4. Schematic of the transition. The ordinate is free energy, and the abscissa is a configu- 
rational coordinate. Stronger excitation causes both a decrease in the relative energy difference 
between the ICM and CM phases and a decrease in the height of the energy barrier. At high 
excitation the change in configurational coordinate is a 1/4 cycle of a low momentum spin wave. 



14 




Supplementary Figure 1: Time-dependence for selected lower excitation 
fluences. (a) The time-resolved changes in the CM and ICM peaks for 14 
mJ/cm 2 . The arrow indicates the time where the relative intensity starts to 
change. Inset: long time behavior, (b) The ratio of the ICM to the CM 
intensity as a function of time for 7 mJ/cm 2 and 14 mJ/cm 2 excitation levels. 
The solid lines are fits to the delayed bi-exponential function described in the 
text. 
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Supplementary Figure 2: Fluence dependence of the ICM/CM ratio at 100 
ps after excitation. Note that for a complete transition over the entire probed 
volume the ICM/CM ratio would be infinite. 
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Supplementary Figure 3: Time-dependence for the commensurate peak 
intensity at a temperature of 150 K and a fluence of 22 mJ/cm 2 . 



16 



